14.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

  • Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  • Independent variables: Study, GA at sampling
  • Model fitted: LME or GEE model
          log(hormones) ~ Study + GA at sampling + Study:GA at sampling
          LME only: Random intercepts for Study ID and Plate
          GEE only: Independence correlation structure

14.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 21.782244 1 <0.0001 ***
Study 1.759986 1 0.2 ns
GA_samp 37.833421 1 <0.0001 ***
Study:GA_samp 1.825160 1 0.2 ns
hPL (μg/ml)
Chisq Df Pr(>Chisq)
(Intercept) 107.37637 1 <0.0001 ***
Study 63.37312 1 <0.0001 ***
GA_samp 115.88145 1 <0.0001 ***
Study:GA_samp 71.34390 1 <0.0001 ***
GH2 (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 9.974225 1 0.002 **
Study 26.572090 1 <0.0001 ***
GA_samp 57.571753 1 <0.0001 ***
Study:GA_samp 26.560416 1 <0.0001 ***
Insulin (uU/ml)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 3.706 3.706 0.05 ns
3 GA_samp 1 Inf 0.183 0.183 0.7 ns
2 Study:GA_samp 1 Inf 2.665 2.665 0.1 ns
HOMA-IR
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 4.810 4.810 0.03 *
3 GA_samp 1 Inf 0.432 0.432 0.5 ns
2 Study:GA_samp 1 Inf 2.900 2.900 0.09 ns
Glucose (mmol/L)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.138 0.138 0.7 ns
3 GA_samp 1 Inf 0.001 0.001 1.0 ns
2 Study:GA_samp 1 Inf 0.190 0.190 0.7 ns

14.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)
Var Effect P
(Intercept) 2.14 (1.24, 3.03) <0.0001
StudySTOP 0.63 (-0.30, 1.56) 0.18
GA_samp 0.18 (0.12, 0.24) <0.0001
StudySTOP:GA_samp -0.04 (-0.10, 0.02) 0.18
hPL (μg/ml)
Var Effect P
(Intercept) 2.53 (2.04, 3.01) <0.0001
StudySTOP -2.02 (-2.52, -1.52) <0.0001
GA_samp 0.17 (0.14, 0.20) <0.0001
StudySTOP:GA_samp 0.14 (0.11, 0.18) <0.0001
GH2 (ng/ml)
Var Effect P
(Intercept) -1.04 (-1.70, -0.38) 0.002
StudySTOP -1.77 (-2.45, -1.08) <0.0001
GA_samp 0.16 (0.12, 0.20) <0.0001
StudySTOP:GA_samp 0.12 (0.07, 0.16) <0.0001
Insulin (uU/ml)
Var Effect P
(Intercept) 3.45 (2.16, 4.73) <0.0001
StudySTOP -1.31 (-2.64, 0.02) 0.05
GA_samp -0.05 (-0.13, 0.04) 0.27
StudySTOP:GA_samp 0.07 (-0.01, 0.16) 0.10
HOMA-IR
Var Effect P
(Intercept) 2.16 (0.83, 3.50) 0.001
StudySTOP -1.55 (-2.94, -0.16) 0.03
GA_samp -0.06 (-0.14, 0.03) 0.21
StudySTOP:GA_samp 0.08 (-0.01, 0.17) 0.09
Glucose (mmol/L)
Var Effect P
(Intercept) 1.66 (1.30, 2.02) <0.0001
StudySTOP -0.07 (-0.43, 0.29) 0.71
GA_samp 0.00 (-0.02, 0.03) 0.84
StudySTOP:GA_samp -0.01 (-0.03, 0.02) 0.66

14.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

14.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 51.11441 8.906805 1114.07386 36.31267 71.94962
STOP 10 62.65847 4.432004 69.44323 54.41322 72.15314
SCOPE 16 150.41614 10.759690 72.33237 130.42768 173.46791
STOP 16 142.93995 12.259227 148.33938 120.65629 169.33912
hPL (μg/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 67.92125 6.0306146 1668.81645 57.06556 80.84205
STOP 10 37.36817 0.9095297 93.09186 35.60500 39.21865
SCOPE 16 187.35085 4.6396951 97.63587 178.36571 196.78863
STOP 16 242.41586 8.6131317 388.30674 226.05960 259.95557
GH2 (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 1.7651770 0.2079050 779.4183 1.4008017 2.224333
STOP 10 0.9757305 0.0225094 193.2049 0.9323297 1.021152
SCOPE 16 4.6213763 0.1093036 180.9841 4.4106582 4.842161
STOP 16 5.1643805 0.2209490 1235.9833 4.7485975 5.616569
Insulin (uU/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 19.77659 4.6279509 Inf 12.50145 31.28545
STOP 10 11.15166 0.3842101 Inf 10.42348 11.93070
SCOPE 16 14.96218 0.6124863 Inf 13.80862 16.21210
STOP 16 13.12747 0.9877112 Inf 11.32756 15.21337
HOMA-IR
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.004074 1.2160703 Inf 3.107908 8.057111
STOP 10 2.360224 0.0858766 Inf 2.197770 2.534686
SCOPE 16 3.587401 0.1563751 Inf 3.293639 3.907364
STOP 16 2.735166 0.2179779 Inf 2.339631 3.197568
Glucose (mmol/L)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.395297 0.3344260 Inf 4.778085 6.092238
STOP 10 4.773670 0.0366162 Inf 4.702440 4.845979
SCOPE 16 5.475695 0.0719187 Inf 5.336536 5.618483
STOP 16 4.691147 0.0561450 Inf 4.582385 4.802490

14.2.5 Comparisons

Hormone Contrast GA_samp Ratio of Geometric means (95% CI) P-value
PRL (ng/ml) STOP / SCOPE 10 1.23 (0.89, 1.69) 0.2
PRL (ng/ml) STOP / SCOPE 16 0.95 (0.84, 1.07) 0.4
hPL (μg/ml) STOP / SCOPE 10 0.55 (0.46, 0.65) <0.0001
hPL (μg/ml) STOP / SCOPE 16 1.29 (1.21, 1.38) <0.0001
GH2 (ng/ml) STOP / SCOPE 10 0.55 (0.44, 0.70) <0.0001
GH2 (ng/ml) STOP / SCOPE 16 1.12 (1.02, 1.22) 0.01
Insulin (uU/ml) STOP / SCOPE 10 0.56 (0.35, 0.90) 0.02
Insulin (uU/ml) STOP / SCOPE 16 0.88 (0.74, 1.04) 0.1
HOMA-IR STOP / SCOPE 10 0.47 (0.29, 0.76) 0.002
HOMA-IR STOP / SCOPE 16 0.76 (0.64, 0.91) 0.003
Glucose (mmol/L) STOP / SCOPE 10 0.88 (0.78, 1.00) 0.05
Glucose (mmol/L) STOP / SCOPE 16 0.86 (0.83, 0.89) <0.0001

14.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

14.2.7 Plots for GA